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Abstract. In this work, we solve a long-standing open problem: Is it true 
that the convergence rate of the Lions' Robin-Robin nonoverlapping domain 
decomposition(DD) method can be constant, independent of the mesh size hi 
We closed this twenty-year old problem with a positive answer. Our theory is 
also verified by numerical tests. 
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1. Introduction 

Domain decomposition (DD) methods are important tools for solving partial dif- 
ferential equations, especially by parallel computers. In this paper, we shall study 
a class of nonoverlapping DD method, which is based on using Robin-Robin bound- 
ary conditions as transmission conditions on the subdomain interface. The idea of 
employing Robin-Robin coupling conditions in DD methods was first proposed by 
P.L. Lions in |24) . In the past twenty years, there are many works on the analysis 
and applications of this DD method: Despres [5] , Douglas and Huang [TH [T3] , 
Deng [ilT], Du [M], Gander et al [20l [21], Guo and Hon [22], Discacciati [9], 
Flauraud and Nataf [16], Gander [HI [H] , Qin and Xu [26l [23 [28] , Discacciati et 
al. |1U] , Lui [2S] , and Chen et al. [H [3] . We should say that the list is far from 
being complete. 

By comparison with other DD methods. Lions' DD method has several advan- 
tages. The iterative procedure is simple and much more highly parallel than others. 
Because it employs Robin conditions, the method is specially suitable for solving 
Helmholtz and time-harmonic Maxwell equations. There exists a lot of works in 
this direction, cf. [i [H [H HJ for details. 

Lions' Robin- Robin DD method was proposed in 1990 [25, see Definition 11.11 
below (without Step 5). The convergence (without any rate) is shown in [24[ [29] . 
Later, the convergence was improved to a geometric convergence [131 [HI [H], i.e, a 
rate of 1 — 0{h). It was first pointed out by Gander, Halpern and Nataf in [3D] that 
the optimal choice of relaxation parameter is 7 = 0(ft.^^/^) and the convergence 
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rate 1 — 0(V7i) could be achieved. Recently, Xu and Qin [3^ give a rigorous analysis 
on this result and shows that the rate is asymptotically sharp. However, without 
enough knowledge on the method, the two parameters 71 and 72 in Lions' DD 
method are set equal, see Definition 11.11 below, by researchers in above references. 
Thus, the rate of l — 0{\/li) is generally believed optimal for the Lions' DD method. 

This paper answers this long-standing open problem: Is it possible to achieve 
a rate of 1 — C for some some constant C > independent of the mesh size hi 
We give a positive answer. Yes, the constant rate of convergence is achieved by 
well-choosing three parameters in the Robin- Robin DD method, 71, 72 and 9, in 
Definition ll.il Roughly speaking, the optimal choices are 

2f — 1 

71 = 0(1), 72 = 0{h-'), and = 

where i w 1 is the ratio of spectral radii of two Dirichlet-Neumann operators on 
two subdomains. It is shown in this paper, by three types of analysis, that the error 
reduction rate of the DD method is optimal, 1 — C. 

Next, we introduce the Robin- Robin DD method through a simple model prob- 
lem. We solve the following model problem in 2D, which is decomposed into two 



subproblems (cf. Figure [T|): 
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where F is an interface separating ili and 512, and n is an outward normal vector of 
r^i at F. The DD method can be applied to general elliptic PDEs, general domains 
and multiple subdomains, cf. [51 129]. 




Figure 1. A domain is decomposed into two subdomains. 
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The Dirichlet and Neumann interface conditions on F in (|l.ip are combined into 
two Robin interface conditions: 

du dw 

(1.2) ^iu+ — ^ ^iw + — = gi on r, 

on on 

du dw 

(1.3) 72" - 7^ = 72W - = 52 on r. 

on on 

Here we allow 71,72 to be any positive constants. For example, when 71 is ar- 
bitrarily close to zero and 72 is close to infinity (but the linear systems would 
become near singular), the method would be reduced to the Dirichlet-Neumann 
DD method. The past researchers all set 71 = 72 = 7 in the Robin interface condi- 
tions, i.e., the two parameters are simultaneously large or small. By selecting two 
parameters correctly, using the original Lions' DD method, this Robin-Robin do- 
main decomposition method should be better than all existing Dirichlet-Neumann, 
Neumann-Neumann and Robin-Robin domain decomposition methods. 

Let Vi = iJQ(i7)|o;. Later, Vi also denotes the restriction of the finite element 
space of grid size h on the two subdomains fli . By (11.21) , we do an integration by 
parts on 57 1 to get 

[ givds — [ (— \- jiu)vds — [ (Vu • Vw -|- Aww)(ix + 71 [ uvds 

Jr Jr on J^^^ 7r 

= / (Vu • Vu — /v)(ix + 71 / 
jQi Jr 



uvds. 
T 



Thus 



where 



ai{u,v) + 'yi{u,v) = (/, t;)oi + {91, v) Vv eVi, 

ai{u,v) = / Vu • Vudx, i = l,2, 
{f,v)n, = / fvdyi, i = 1,2, 
{u, v) — uvds. 



Similarly, by (jl.3p and an integration by parts on fl2, it follows (noting that n is 
an inward normal vector to 5^2) that 



wvds 



I g2vds — I (72W " ^^)vds — ( {Vw ■ Vv + Awv)dx + J2 [ 
Jr Jr on J^^ Jr 

w = / (Vw • Vw — fv)dx + 72 / wvds. 
Jn2 Jr 

This way, we get the second variational problem on D,2- 

a2{w,v) +-i2{w,v) = (/,w)o2 + (32, "> Vu e V2. 

Definition 1.1. (The Robin-Robin DD method.) Given gi{= 0) on T, a se- 
rial version domain decomposition iteration consists the following five steps (m = 

0,1,...;.- 
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(1) Solve on fli for u™; 

(1.4) ai{u"\v)+j,{u"',v)^{f,v)n, + {gT,v) yveV,. 

(2) Update the interface condition on T: 

(1.5) 52"--5r+(72+7i)"". 

(3) Solve on ^2 for w"^ : 

(1.6) a2{w'^,v)+j2{w'^,v}^{f,v)n, + {gT,v} Vi; e ^2- 

(4) Update the other interface condition on T: 

(1-7) 9T^-9T + {li+l2)w"^. 

(5) Get the next iterate by a relaxation: 

(1.8) gT+'^0gT + il-9)gT. 

The rest of paper is organized as foUows. In section 2, we shall show that 
although the Robin-Robin DD method cannot achieve the geometrical convergence 
rate at the continuous PDE level, but it does at the discrete level. In section 3, 
we shall give an explicit convergence rate of the DD method on uniform meshes. 
In section 4, we shall extend our method to more general quasi- uniform meshes. 
Using the Dirichlet-to-Neumann operator, we shall prove that the Robin-Robin 
DD method is optimal. Finally, in the last section, we shall present some numerical 
results to support our theory. It is seen from our numerical implementation that 
this DD method is better than Dirichlet-Neumann DD method and one-parameter 
Robin-Robin DD method. 

2. A Von Neumann analysis 

In this section, through a simple model problem, we shall show that for the new 
DD method it is not possible to get the geometrical convergence rate strictly less 
than one at the continuous level, but it is possible at the discrete level. 

Let us assume that fli = [— tt, 0] x [0, tt] and Q2 = [0, tt] x [0, tt], and it is enough 
for us to assume that / = so that the true solutions of Equation (jl.ip vanishes. 
Now if gi = gi sinfcy on F, from Equations (jl.ip and (|1.2p . the solution on fti is 

M = 'usinh(fc(a; -f- 1)) sinfc?/, where u— 



71 sinh k + k cosh k 

If 92 = .92 sin ky on F, from Equations (jl.ip and (jl.3l) . the solution on D,2 is 

g2 

w = wsinh(fc(x — 1)) sinfcy, where w — — ; -— . 

72 smh fc -I- fc cosh fc 

By Definition ll.il if the initial error is = ffi" sinfcy on F, then 

= sin ky, where q^ = g?" ( ^\ , - 1 

^' V7i + fccothfc 
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Then, by and ([TJl) . 

72 + 71 



~gT = ~gT sin ky, where 5^-52™ 
Finally, after the relaxation step (|1.8p . 



72 + fc coth fc 

g™+i^5™+isinA;2/, 



where 

and the factor 

(2.1) p^9+il-9)( 

^ ^ ^ ^V72 + ^cothA: yV7i + fccothfc 

Now for the fixed parameters O<0<1, 7i>0 and 72 > 0, if fc tends to infinity, 
then 

(2.2) p^l-2(l-g) ^'|^' ^0. 

k 

Therefore in the continuous case, it is impossible to get the convergence rate in- 
dependent of the frequency (or the wave number) k. On the other hand, if k is 
bounded by 1 < fc < K, we may obtain the convergence rate p , which is indepen- 
dent of k (but dependent on K) , through choosing the three parameters 71 , 72 and 



Lemma 2.1. If a and b are two non-negative constants, then the function 
(2.3) p{9) = max{|6l - (1 - e)a\, \9 - {1 - 9)b\} 

attains the minimum value at 9n = .-, . . 




Figure 2. Graphs of \9 ~ {1 - 9)a\ and \9 - {1 - 9)b\, b> a>0. 



Proof. Without loss of generality, we assume b > a. Both terms in (|2.3|) are piece- 
wise linear functions. We plot them in Figure [5J The minimal value is attained at 
the point A, where two lines intersect: 

9-{l-9)a + 9-{l-9)b^0. 
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That is 9 = 00= 2^^, and 



\9o - (1 - 0a)a\ = \eo - [1 - 9o)b\ 



a + b 



(2.5) maxwfz) = -^r, where rj ■ 



2 + a + b 

So we get the lemma. □ 

Lemma 2.2. For any z > 0, the function 

(2.4) uj{z) = — — 

72+2 z + 7i 

attains the maximum value at zq = ^7172; 

fcll! where 
Proof. The derivative of a;(z) is 

/ORN V A 2(71 +72)(7l72 - z^) 

(2.6) UJ (z) = — -r- -7^. 

(^2 + 71)^(^ + 72)2 

So uj(z) monotonically increases when z < zq and monotonically decreases when 
z > Zq. In particular, then the minimum value of uj{z) on an interval [zi,Z2] is 
attained at one of the end points: 

(2.7) min Cli(2;) = min{a;(zi), ^(2:2)}. 

By (j2.6l) . uj{z) attains the only global maximum value at zq — ^7172: 
, . 72 - V7172 V7172 - 71 iv- 1)^ 

■^(-^O) = = -7r. 

V7i72 + 72 V7i72 + 71 (jy + i) 
The lemma is proved. □ 

From Equation (|2.ip . 

X 72 — fc coth fc fc coth /c — 71 , , , 

p = 6* - 1 - e/U^ ^ — -4 = 9-(l- 9)oj(kcothk). 

' ^ 72 + fc coth k 71 + fc coth fc \ I \ I 

If 7i and 72 are chosen such that 71 < coth 1 and 72 > cothiC, then cj(fc coth fc) > 

0. By Lemma [23 

IpI <max{|0|,|0-(l-0)c.(zo)|}. 
Applying Lemma |2. 11 we may select 

(2.8) ""=2^ - H<W<i 

Remark 2.1. If rj > ^{^\ , we may just set 9 — -^j and \p\ is also less than i. 
Moreover, this bound can be improved further if we carefully estimate the minimum 
value ofuj(z). 

Remark 2.2. The constrain 71 < cothl can be relaxed. Actually, if ii > cothl, 
then 

\p\ < max{|0 + (1 - 9)C\ ,\9-{l- 9)oj{zo)\} , 
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where ( = ^7+f§^FT- ^^^^^ set 9 = if w(zo) < C and set 9 = 2+If('zo)-C 

iVote t/iat i^izo) < 1 and 2+lf(°zo"|-C — which is also independent of K . 

The Von Neumann analysis shows that the Robin- Robin DD does have a constant 
rate of convergence, independent of the frequency number k or K. But the selection 
of the two parameters depends on K. The limit case indicates that the method 
deteriorates to, i.e., 72 ~ cxd, a Robin-Dirichlet DD method. 




3. Convergence on uniform grids 

In this section, we analyze the Robin-Robin DD method on uniform grids. In 
this case, we give explicit eigenvalues of the iterative matrix, and show the optimal 
rate of convergence. 

We post a uniform grid of size h — l/(2n) on the domain ft = [0, 1]^, the unit 
square, shown in Figure ID Then, we subdivide the domain into two, as shown 
in Figure [3l We give two numberings of nodal values of the C^-Pi finite element 
functions. One numbering is on the interface F. The other one is within each 
subdomain, ili and r22- When numbering the nodes in 172, we go from right to left 
so that the nodal index is symmetric to that on the left domain fli. 



2n-l 





X2n-1 




X3 , , 




X2 ^ ^ 




Xi ^ ^ 




^2 



(2n - l)n(last) 



1 


/ 


2 


2 


1 


1 




^2 



2n- 1 



Figure 3. Nodal basis numberings, on F and on fii n il2- 



Let Mr and Ar be two tridiagonal (2?i — 1) x (2?i — 1) matrices: 

/4 1 \ /4 -1 \ 



Mr = 



1 4 



At 



1 4j 



-1 4 



-1 4 J 
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Here Mr is just the mass matrix of the imier product (•, •). Let Rh be the (2n — 
1) X (2n — l)n matrix representing a restriction operator on F: 

(3.1) Rh = (02n-l, • • • , 02n~l, hn-l)- 

The stiffness matrix of the bihnear form ai(-, •), under nodal basis, (and a2(-, •) too) 



IS 



Ah 



_ Ao-RhArRh, 

where the matrix Aq is the stiffness matrix of size (2?! — 1)?!, for the Laplace 
operator on a (2n) x {n + I) uniform grid with zero Dirichlet boundary condition. 
^0 is same as the matrix of standard five-point finite difference matrix, which has 
the eigen-decomposition O [23] : 

(3.2) Aq = ($„ (g) $2n-l)'^(A„ (g) hn-l + /„ » A2„_l)($„ (g) ^2n-l), 

where Am denotes an diagonal matrix whose (i, z)-th entry is 



(m) 



4 sin 



(3-3) - 2(^+1) 

and $m denotes an orthogonal matrix defined by 



/sin^\ 

m+1 \ 



(3.4) 



with 



A-n 



m+1 

2i7T 

m+1 



ysin ^ / 

^ m+1 / 

Bkxk is defined with 



Here in p.2|) . a tensor product matrix Cmkxmk 
the (i, j)-th entry 

Cij = Aii jiBii/ ji/, where i = [i' — l)k + i". 

In Definition 1 for (|1.4p . the error e™ = u — u'" satisfies the equation: 

ai(C,«) + 7i(C,^') = (e^,«) V«eli. 
Here e™ = gi — <7i" is the error. In the matrix-vector form, 

= (A + ^^RlMrRhT^RlMrE^. 
Here i?™ is the vector representation of e™. Therefore, by l|1.6p . 
(3.5) = (-/ + (72 + 7i)i?;»(^h + iiRlMTRhr^RlMv) E^^. 

Symmetrically, by and pT^ . 



(3.6) 



Here, for simplicity, we denote the error reduction matrix by 
(3.7) C^, - (-/ + (72 + -/i)Rh{Ah + -ikRlMvRhr^RlMv) , 
Finally, by (|1.8p . one Robin- Robin DD iteration reduces the initial error i?'" to 



1,2. 



(3.8) 



Tpm+l 
91 



[ei + {i-e)c^,c^,]E";^. 
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We find the eigenvalue range of this error reduction matrix, via common eigenvec- 
tors of all matrices. 

Lemma 3.1. The error reduction matrix p.8|) can be diagonalized by $2n-i defined 
in p.4p . That is, 

(3.9) $2„-i[e/ + (1 - 0)C^.C^,]<fl,^, = diag(0 + (1 - e)cj), 
where in the j-th diagonal element, 

(3.10) c, = ^'"^-^^ ■ 12^1^. 

"fia-j + bj 72%- + bj 

Here in (|XTU)) . 



h . (2n-l) 



(3.11) flj = XMrJ^j, ^Mr,j — h — —X- 

(3.12) b, = l-XA,/X,, ^Ar,j = l + \xf'-^\ 
where A^-^" is defined in (|3.3p and 

(3.13) A, = ^ Vsin^^(Ar) +Af -^V^ 

z— 1 

Proof. In p.5|) . by the Sherman-Morrison- Woodbury formula, 

= (Ao + RI{-At + 7iMr)i?,,)-i 

= - A-ii?^((-^r + IxM^r^ + i?;,A-ii?;^)-ii?„A-i. 

Now letting Bq = BhA^^ R-^ , we have 

RhiAh + -fiRlMvRhr^Rl = Bo- BoH-Ar + IiMtY^ + B^y^B^. 

By (|3.ip and (|3.2p . notice that ($„ (g) <i>2„-i)i?^ = ^l"' <i>2ri-i, we can compute 
Bo: 

Bo = (0l") ® $2n-^l)'^(A„ (g) /2„-l + /„ «) A2„_l)(0l") ® $2n-l) 
n 

- + A2„-l)-l$2„-l 

= <fL-l (^E(0i;))^(Af"-^V2„-l + A2„-l)-^) $2„-l 
= *L-lAo*2„-l, 

where is the i-th entry of vector (j^n^ defined in (13. 4p . and Aq is a diagonal 
matrix, whose (j, j)-th entry is define in p.l3p . The matrices on V are diagonalized 
as Mr = ^'L-i diag(AAfr j)"J'2n-i and Ay = diag(AAr,j)*2n-i, where Ajv/r j 
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and XAr,j are defined in (|3.1ip and p.l2|) . respectively. Thus combining last two 
equalities, we get 

RhiAh+iiRlMrRhr^Rl 
= Bo -Bo($^„_i(-diag(AArj) +71 diag(AMr j))"^*2n-i + So)"^So 



2n-l 
T 

2n-l 



Ao-A^{(-diag(AAr,j)+7idiag(AMrj)) ^ + Aq} ^ 



Aq ^ - diag(AAr j) + 71 diag(AA/r,j) 

By en, 

C71 = + (72 + 71) [(A;^' - diag(A^^ j)) diag(A^^^_^.) + 71/)] ')$2n-i 

/ — 1 + 72AMr j-^j + '^Ar,jXj 



= 3'2„-idiag 



1 + ll^MrjXj ^ ^ArjXj 
$2«-l- 



2n-l 



-C-idiag(^^ 
In the same fashion, it follows that 



C72 = *2«-i diag 



/ 7iaj - bj 



2n-l- 



Thus (|331) follows. □ 

In next lemma, we estimate the eigenvalue Cj in the reduction matrix, p.lOp . 
Lemma 3.2. (3aj — hj) is monotonically decreasing, i.e., j = 1, . . . , 2n — 2, 



(3.14) 



3a, — bj > 3a,+i 



Oj + l. 



Proof. We rewrite the Xj (in and bj) in a symmetric form so that each i-term is 
a decreasing function of j (the original term is not.) 

2 sin2(i7r/(n + 1)) 



1 2 ^ 



n + 1 2 ^ 4sin2(j7r/(4n)) +4sin2(i7r/(2n + 2)) 
sin^((n + 1 - i)7r/(n + 1)) 



+ 



4sin2(j7r/(4n)) +4sin2((n + 1 - i)7r/(2n + 2)) 
1 A Sin2(i7r/(n + l))(2sin2(j7r/(4n)) + l) 



n + 1 ^ 4sin^(j7r/(4n)) + 4 sin^(j7r/(4n)) + sin^(i7r/(n + 1)) 
To shorten expression, we introduce two more notations 



e • 2 JT^ 
^] sm — , 
4n 



(3.15) 

(3.16) 6lj = (1 + 2Q{1 + 3h + - 2h^j). 

By (pHI) and (151^ . we have 



(3.17) 



3aj - b^ 



sin^(i7r/(n + l))6'j 
+ 1 ^ 4e|+4e,+sin2(iV(n + l))' 
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We show that each term is a decreasing function of ^j. That is, each term 

^ (2e + l)(2(l + M-(l-fe)(l-20) 
•^'^^ 4e2+4e + sin2(^7r/(n + l)) 

is a decrease function of ^, for ^ e (0, 1). By the quotient rule, 

(4(1 + h) + 8(1 - h)Oi'ie + 4e + sin2(«7r/(n + 1))) 



mo 



(4^2 + 4^ + sin2(i7r/(n + 1)))2 
((1 + 3h) + 4(1 + h)^ + 4(1 - h)e)i8i + 4) 



(4^2 + 4^ + sm^{iTr/{n + 1)))2 

The combined numerator is 



4(1 + h) cos^ + 8hj - 1^8(1 - /i) cos2 — -j-^ + 16/ij ^ - {32h)C < 0. 

As each term fi{^j) is desecrating with respect to j, the sum is a desecrating 
function of j. We prove the lemma. □ 

We will find a bound for the biggest term (3ai — 5i), among all (Soj — bj), in 
order to bound the Cj in p.lOp . One can prove that, for all n > 1, 

(3.18) 3ai-bi < -Ih^/ie, 

where oi and hi are defined in (|3.1ip and p.l2p , respectively, and h — 1/ (2n) . But 
our proof for (|3.18p is lengthy and tedious. In this paper, we prove a worse bound 
only, in the next lemma. 

Lemma 3.3. If n > 11, then (cf. (jXTS)) } 

(3.19) 3ai - 6i < ~0.049/i < -7h^/16, 
where ai and hi are defined in p. lip and p.l2p . respectively. 

Proof. By ((3T7l) . with the notations defined in (|31T|) . ((3T2t and (|3T3)) . 

(3.20) 3ai - &i + 1 = (1 + 3/i + i^A^ ""'^)Ai. 
We estimate an upper bound for 



Ai = 



2 eos2^sin^ji, 

+ 2(n+l) n + l^sin^^^+sin^^ 

n 2sin2^A/ 1 + \ 



n+1 n + 1 Z_.^3in2_^_^3j^2^ 
l-2nsin2JL 2(l + sin2^) 



3.21 =1 4n__^ 



sm 



2 ji_ 

4n 



n+1 n+1 ^ sin^ tj^^ + sin^ ^ 

1=1 2{n+l) 4n 
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As (sinx/x) is a decreasing function of x on (0, 7r/2), we have 



V ^ > V > V- 



11 1 
> Y 7^ > 0.33462. 

z— 1 

Substituting the estimate into the expression of Ai , 

^ 1 + 2(1 + 2 sin^:!^)- 0.33462 - 2nsin2^ 1.55726 

Al < 1 ; < 1 — . 

n + 1 71+1 

By dSjO]), if n > 11, 

. . 3 7r2 ^ , 1.55726, 

< 1 - 0.049/1 < 1 - OMh^ < 1 - 7/iVl6. 
We proved the lemma. □ 

With the expUcit eigenvalues of the reduction matrix and their bounds, we can 
easily choose a set of parameters 71, 72 and 9, to get a constant rate of reduction, 
independent of mesh size h. 

Theorem 3.1. Let 71 — 1, 72 = 64/i^i and 9 ~ 3/7 in Definition ] 1.1[ The error 
reduction factor (for the Pi finite element on uniform grids shown in Figure ^ is 
bounded by 1/7, independent of the grid size h, 

II „m+l II ^ ^ II „m II 

11^91 lk=^(r) < yl|egJ|L2(r). 

Proof We will apply Lemma O By (pUj) and (031), > 0. By (|XTi)) . (|XT5)) 
and (P+a . 

3aj ~ bj < 3ai — 61 < — 7ft-^/16, 
(3.22) bj > 3aj + 7h^/m > 0. 

By (Eini), 

1 — fej/flj 64/i^i — bj /oj 
'^^ ~ 1 + bj/ttj 64/1-1 _,_ fj./fj^. ■ 

We let z = bj/uj > in Lemma [521 The critical point is (may be outside the bj/aj 
range) 

We find the two end points of possible z. First, by p.l3p . 
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Thus, by ([3lT|) . ([332]) and (fS^]) . 

< (/i - ^ • 0) • 1 = /i, 

.,<i-a + i.o)^l.I. 

In the first inequality, we used p.2ip tliat < 1. Wc find one end point for z: 

aj - h/12 2h 

For tlie other end point, by (I3.22p . 

bi 7h^/m 7/iVl6 7h 

- > 3 + — > 3 + ^3+ — =zi. 

aj aj n lb 

By Lemma 12.21 the range of Cj is between its values at z — zi, zq, Zr. We note that 
zi < zq < Zr here. At each point, we need to apply [2?2] again for h varying. But we 
can find some rough (but good enough) bounds at each point, directly. 

2184975 1 

At z = zi: -0.50098... = < c,- < — = -0.5. 

4361329 - ^ 2 

32- 129%/2 

At z = Zq: -1 < Co < V = -0.7015... 

32- 129\/2 

Hence the value of Cj is always strictly between —1 and —1/2. When 9 — 3/7, 
we get, 

(3.23) + (i_0)c^. > 3 + 4^_l) = _l, 



(3.24) + (I - 0)cj < - + -{--) 



3 4 11 
7^7^ 2' ~ 7' 



This gives the error reduction factor. □ 

By (|3.23p and (|3.24p . we can get the following result for a general relaxation 
parameter 9. 

Corollary 3.1. Let 71 = 1 and 72 = 64h~^ in Definition ] The error reduction 
factor p Jar the Pi finite element on uniform grids ( shown in Figure is 

jl-29, 0<9<3/7, 
^~ [(361- l)/2, 3/7<0<l. 

That ts, |le^+i|lL^(r) <p||e^|U2(p). 
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4. Convergence on general grids 

In this section, we consider the convergence behavior of the Robin-Robin DD 
method on general quasi- uniform meshes. By the algorithm in Definition 11.11 for 

i = 1,2, 

a,{eT,v)+^,{eT,v)^{e^\v), 

where the errors are defined by 

£™ = 5^-c,™, e'i" = u-u", and ^ w - w"^ . 

Let 5*1 and S2 be the standard Dirichlet-to- Neumann operators, cf. [291I3Q]. The 
error e" {i = 1, 2), restricted to the interface F, satisfies the relation 

(4.1) er = (7. + 5,)erir. 
Using the first interface update (|1.5|) . we have 

(4.2) e2" = -er + (7i+72)erir. 
For the second one, by (|1.7p and (|1.8I) . 

- 6eT + (1 - 6)[-e^ + (71 + 72)6^" |r] 
= fer + (1 - ^)h(72 + S2)e^\T + (71 + 72)e^|r] 
= fer + (l-0)(7i-^2)erir. 
By dm]), (gj]), we have 

6^+' = 6er + (1 - ^^)(7i - ^2)(72 + S2)-'e^ 

- + (1 - ^^)(7i - 52)(72 + S2)-\i2 " 5i)er|r 
= [e + {l- ^?)(7l - ^2)(72 + ^2)^^(72 - ^i)(7i + Si)-^]eT. 
Let us represent the iteration by 

where 

(4.3) R = e-{1- e)T, 

(4.4) T = (^2 ~ 7i)(72 + 52)-i(72 - 5i)(7i + Si)-\ 
Next, we give an convergence analysis for this DD operator R. 

4.1. Symmetric case: Si — S2{'-= S). Let z be an eigenvector of the symmetric 
operator 5* (cf. [521130]) corresponding to the eigenvalue A^. By (|4.4p . z is also an 
eigenvector of the symmetric operator T. 

^ [6 - {1 - e)Lo{X,)]z. 

It is known [30] that 

e [co,Coh-^]. 
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Now if we choose 

(4.5) 0<7i<co, and "f2>Coh'\ 

by Lemma [221 we get 

Then we bound the spectrum of the symmetric operator i?, (t{R), as 

^(i?)c[.-(i-.)|^,.]c[4iL 

when choosing the parameter 9 = 1/3, cf. Remark 12.11 That is, the convergence 
rate is bounded by 1/3, independent of the mesh size h, when choosing parameters 
by C31) . 



4.2. Nonsymmetric case: Si « S'2. In this case, there exist two positive constant 
< s < 1 and t > 1, independent of the grid size h, such that for all w G V^jr (cf. 
for details): 

(Al) s{Siv, v) < {S2V, v) < t{Siv, v). 

Si{i = 1,2) are symmetric and positive definite(SPD). Let Aj be the minimum 
eigenvalue, and Xi the maximum eigenvalue of Si. In this subsection, we assume 
that the parameters are chosen to satisfy 

(A2) < 7i < min{Aj, Aj}, and 72 > 3max{Ai, A2}. 

The parameter selection is similar to that in the symmetric case, (j4.5l) . 

Lemma 4.1. The condition (lAlll has another version 



(4.6) \{S^'v,v) < {S^\,v) < ks-\,v). 
Proof. Replacing v hy S^ in (IA1|) . 

s{v,v) < {s;h2s;K,v) < tiv,v). 

This inequality implies that the spectrum of the SPD operator Si '^8281 ^ is within 

~ 1 ~ 

So the spectrum of its inverse, S2 Si is inside [t~^^s~^]^ i.e., 

j{v,v)<{S^S^'S!v,v)<kv,v). 

(I4.6P follows after replacing v by S~^v. □ 

To find the spectrum of DD operator T in (|4.4p . we introduce a symmetric 
operator 

(4.7) f = (71 + 5i)-5(72 - Si)HS2 - 7i)(72 + S2)-Hj2 - 5i)^(7i + ^i)"^. 
This operator is similar to the nonsymmetric operator T, defined in (j4.4[) . 

Lemma 4.2. If Assumption (jA2[) is satisfied, then T is SPD. 
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Proof. T is symmetric because 

= f. 

Notice that (S'2-7i)(72 + 52)"^ = /-(7i+72)(72 + 5'2)"^ Its minimum eigenvalue 
is 

l-(7i+72)(72 + A2)-i-^^, 

72 + A2 

which is positive by Assumption ()A2|) . Similarly, the minimum eigenvalue of (72 — 
'S'i)(7i + Si)~^ is (72 — Ai)/(7i + Ai), which is also positive by Assumption (jA2p . 
Now for any v S Vi\r, we have, denoting v — (72 — Si)i{'^i + Si)^^v, 

{Tv, v) = ((S'2 - 7i)(72 + S2y^v, v) 

. A2-71 
72 +A2 

= ^^^((72-5i)(7i+5i)-V«) 
72 +A2 

A2 - 71 72 - Ai 

= — • —[v,v). 

72 + A2 71 + Ai 

~ • Ao — 71 72 ~ Ai 

It means that the minimum eigenvalue of T is greater than — > 0. 

72 + A2 71 + Ai 

That is to say, the symmetric operator T is also positive definite. □ 

We find an upper bound of the spectrum of SPD operator T next. To this end, 
we rewrite T as 

(4.8) f^flf,f2, 
where 

(4.9) fi = (^2 - 7l)(72 + S2)-Hl2 - S2){jl + S2)-\ 

(4.10) f2 = (71 + ^2)^(72 - ^2)"^ (71 + ^i)~^(72 -Si)i. 
Lemma 4.3. // (ET|) and ([M]) hold, then, for the t defined in IXTjl . 

(4.11) ((72-^2)"'(7i + ^2)t;,«) < {2t-l)((j2-Si)-\'ji+Si)v,v). 
Proof By jJT]) and 

((72 - <S'i)(7i + Siy^v^v) = ((71 +72)(7i + Siy^v,v) - {v,v) 

< t ((71 + 72)(7i + S2y^v, v) - (w, v) 
= t ((72 - ^2)(7i + S2)-'v, v) + {t- l)iv,v). 
We bound the second term next. By the assumption (|A2[) . 

((72 - 5'2)(7i + S2r^v,v) > — — ^{v,v) > ^{v,v) = {v,v). 

71 + A2 ZA2 
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Combining above two inequalities, 

(4.12) ((72 - 5i)(7i + S,)-'v, v) < {2t - 1) ((72 - S2){li + S^T^v, v). 

Applying Lemma |4TT| replacing 5*2 there by (72 — >S'i)(7i +51)^^ and 5*1 by (72 — 
^2)(7i + ^2)-\by 

^^3Y((72 - <S'2)"^(7i + S2)v,v) < ((72 - 5'i)"^(7i + Si)v,v). 

()4.1ip is proved. □ 

Lemma 4.4. // assumptions (jAl[) and (jA2[) hold, then the spectrum of the SPD 
operator T is bounded by 

(4.13) a(T') C (0,2t- 1]. 
Proof, fi is SPD, cf. (gH). The eigenvalues of f 1 are 



^2,j - 71 72 - ^2,j 
72 + A2 J- 71 + ^2,j ' 

where {^2.3} are all eigenvalues of 82- By (|A2[) . 



A, > ^'^ , -"'^ > 0, 

72 + A2,j 71 + A2,j 

A, < ^ < 1. 

71 + A2j- 72 + A2,j 

Then, by and gH]), 

< {Tv, v) < if2v, f2v) = ((71 + 52)(72 " S2)-^i, i) 

< (2t - l)((7i + ^i)(72 - ^i)-'^,^) = (2t - l){v, v), 

where 

W = (71 + S'l)"5(-y2 - Sl)h. 

The proof is completed. □ 



Theorem 4.1. If the assumptions (jAip and (jA2p /loZd, t/ien the spectrum of DD 
reduction operator R, defined in (j4.3p . is bounded, independent of the grid size h: 

r 2t-l 2t-l, 

when 9 is selected by 

91 - 1 

(4.15) = 



2t+l 

Proof By (|4l^ and (|43t . 

a(i?)c [0-(l-0)(2t-l),0], 

where t is defined in (|A1[) . independent of h. Similar to the idea in Lemma 12. 1[ 
(|4.14l) follows after we choose the optimal 6 by (|4.15p . □ 
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5. Numerical test 

For numerical test, we solve the Poison equation (|1.1[) on the unit square [0, 1]. 
The exact solution is chosen 

We choose a; = 1/2 as the domain decomposition interface. We use Pi conforming 
finite element on uniform criss grids, shown in Figured 




Figure 4. A uniform criss grid of size ft, = 1/8. 

First, we do the Robin- Robin iteration (Definition II. ip for problems with dif- 
ferent grid size. The parameters used are 71 = 1, 72 = 64/ ft and 9 = 3/7. The 
iteration stops when Ig™"*"^ ~ 9T\i°° < 10~^^. The number of iteration, the error 
and the order of convergence for the finite element solution are listed in Table 15.11 
We note that there is a superconvergence for the finite element solution in scmi-_ff ^ 
norm. 

Table 5.1. The errors and the iteration numbers, by Definition ll.il 



ft 




ft" 


\ui - Uh\H^ 


ft" 


#DD 


1/4 


0.0027120 




0.203663 




14 


1/12 


0.0000716 


1.65 


0.004456 


1.74 


14 


1/20 


0.0000098 


1.93 


0.000605 


1.95 


14 


1/28 


0.0000026 


1.97 


0.000159 


1.98 


14 


1/36 


0.0000009 


1.99 


0.000058 


1.99 


14 


1/44 


0.0000004 


1.99 


0.000026 


1.99 


14 


1/52 


0.0000002 


1.99 


0.000013 


2.00 


14 



Next, we check our theoretic bounds in Theorem 13. II In (|3.23p and (|3.24|) . if we 
vary Q from 0/7 to 6/7, we can get the following theoretic bounds: 

7 5 3 1 5 8 11 

7' 7' 7' 7' 14' 14' 14' 
We compute the real bounds for these on various meshes, and list them in Table 
15.21 We note that, when Q = 0/7 = 0, the method is reduced to the traditional 
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Robin- Robin DD method (by other researchers, where 71 = 72), which converges 
at a rate of 1 - CVh, cL [SO]- This can be seen in the first column of Table 15.21 



Table 5.2. The reduction rate with different 9 in Definition ll.il 



h\0 





1/7 


2/7 


3/7 


4/7 


5/7 


6/7 


1/4 


0.764 


0.512 


0.260 


0.096 


0.322 


0.548 


0.774 


1/12 


0.865 


0.598 


0.332 


0.115 


0.336 


0.557 


0.779 


1/20 


0.894 


0.624 


0.353 


0.116 


0.337 


0.558 


0.779 


1/28 


0.910 


0.637 


0.364 


0.116 


0.337 


0.558 


0.779 


1/36 


0.920 


0.646 


0.371 


0.116 


0.337 


0.558 


0.779 


1/44 


0.927 


0.652 


0.377 


0.116 


0.337 


0.558 


0.779 


1/72 


0.943 


0.665 


0.388 


0.116 


0.337 


0.558 


0.779 


1/288 


0.971 


0.689 


0.408 


0.126 


0.337 


0.558 


0.779 


1/1152 


0.985 


0.702 


0.418 


0.134 


0.337 


0.558 


0.779 


Corollary 13.11 


1.000 


0.714 


0.428 


0.143 


0.357 


0.571 


0.786 



Finally, we compare the Robin-Robin domain decomposition method with the 
traditional Dirichlet-Neumann domain decomposition method. We code directly 
the Dirichlet-Neumann domain decomposition method, defined as follows. 

Definition 5.1. (The Dirichlet-Neumann domain decomposition method.) 
Given 0) on T, find u"" £ Vi, u'^lr = w"" : 

Find e V2: 

a2{w"-+\v) = if, v)n, - ai(u", v) Vv € V2, 
where v is extended into fli with nodal values. Then 

w"'+^ = Our + {1 - e)w"'+\ 

In Table 3, we list the number of Dirichlet-Neumann domain decomposition 
iterations for the above test problem, for various 9. It seems that no matter how 
to choose 9, the Dirichlet-Neumann domain decomposition method (21 iterations) 
is worse than the new Robin- Robin domain decomposition method (14 iterations). 

Table 5.3. The iteration number for Dirichlet-Neumann DD 
(Definition O) 



h\9 





0.25 


0.35 


0.4 


0.45 


0.5 


0.55 


0.75 


1/4 


88 


24 


22 


25 


29 


33 


38 


78 


1/12 


237 


34 


21 


23 


26 


30 


35 


71 


1/20 


392 


37 


22 


22 


25 


29 


33 


68 


1/28 


548 


38 


23 


21 


24 


28 


32 


66 


1/36 


705 


39 


23 


21 


24 


27 


31 


64 
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